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i We present a comparative analysis of multifractal properties of financial time series built on 

stock indices from developing (WIG) and developed (S&P500) financial markets. It is shown 
how the multifractal image of the market is altered with the change of the length of time series 
f-H , and with the economic situation on the market. We emphasize that the proper adjustment of 

I scaling range for multiscaling power laws is essential to obtain the multifractal image of time 

^ ' series. We analyze in this paper multifractal properties of real financial time series using Holder 

■ /(ct) representation and multifractal-DFA method. It is also investigated how multifractal 

Q^l properties of stocks change with variety of "surgeries" done on the initial real financial time 

series. This way we reveal main phenomena on the market influencing its multifractal dynamics. 
In particular, we focus on examining how multifractal picture of real time series changes when 
i one cuts off extreme events like crashes or rupture points, and how fluctuations around the 

I main trend in time series influence the multifractal behavior of financial series in the long-time 

on ' horizon for both developed and developing markets. 

m 
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1 Introduction 



Multifractality [1-4] is rather a new concept when applied to financial time series. It may be con- 
sidered as the higher order extension of the monofractal analysis used successfully e.g. in analysis 
of persistency level in data - particularly in econophysics. So far it is not very clear what practical 
applications of multifractality in financial time series might be. However, some preliminary directions 
in this topic have already been shown (see e.g. [5]). In monofractal approach, one usually looks for 
the probabilistic (stochastic) behavior of some geometrical or topological properties of signal fluc- 
tuations around the local trend in given time series. One expects power-law relation between the 
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quantity describing such fluctuations (e.g. variance) and the length of time window s along which 
the fluctuation is being measured. A good example of such relation is the one provided in Detrended 
Fluctuation Analysis (DFA) [6-11] where the power law relation between the variance of detrended 
signal F^{s) and the width of time window s reads 

F^{s) - s^^ (1) 

with the Hurst scaling exponent H {0 < H < 1) [12,13]. < H < 1/2 the signal is said to be 
antipersistent, if 1/2 < < 1 the signal is persistent, while the case H — 1/2 corresponds to no 
memory present in signal increments (Brownian motion). 

The relation given by Eq.(l) is usually understood as independence of scaling properties of a system 
from the scale - thus H is constant. However, in many cases it might not be so. There might be a 
finite or infinite number of crossover points sx such that for time scales s < sx the fractal properties 
differ from these for s > Sx- In this case, the fractal structure will be described by the whole set of 
H exponents instead of just one, thus revealing the multiple scaling rules. Using the standard DFA 
procedure one calculates only the leading scaling rule linked to major or more frequent fluctuations. 
Moreover, distortion of the signal from the leading pattern before or below the crossover point sx 
may be very weak, so in order to extract them, one is forced to use more sophisticated technique 
based on the artificial amplification of weak (with respect to the average ones) fluctuations. Such 
philosophy makes the background of the modified DFA proposed by Kantelhardt et. al. and called 
multifractal DFA (MF-DFA) [14]. One replaces in this method the "ordinary" fluctuation function 
F^(s) from Eq.(l) by its q-th moment F^{s) defined as follows 

2N 

F{s,q) = {l/2Nj2[F\s,k)f'y/'' (2) 

k 

for q ^ 0, and 

2N 

F{s,0)=exp{l/4Nj2HF\s,k)]} (3) 

k 

for q = 0, where N counts the number of non-overlapping boxes of size s each for which detrendization 
procedure is performed. 

One obtains in this manner the whole continuous set of Hurst exponents h{q) labeled by different 
qeH. They are related to different level of amplification of small fiuctuation in data. The h{q) 
dependence is decreasing monotonic function of q for stationary signal what basically rcfiects the 
fact that relatively small fluctuations happen more often in this signal than relatively big ones do 
[14]. If so, the so called Holder spectrum or singularity spectrum f{a) [4] of the Holder exponent a 
ca be used where 

a = ^ (4) 
aq 

and 

s{q) = qh{q) - 1 (5) 
while f{a) is determined by Legendre transform [16,17] 

f{a) = q{a)a - s{q{a)) (6) 
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The latter approach makes a useful representation of multifractality. The singularity spectrum in this 
case should have rather regular form of inverse parabolic shape as shown in Fig.2. The width of /(a) 
spectrum measures the multifractality level of the signal. For the pure monofractal signal this width 
converges to one point a — H — h{2) and /(a) = 1, while if the bigger amount of multifractality is 
present in the signal, the width of f{a) spectrum becomes wider. 

We will analyze in this paper multifractal properties of real financial time series using /(a) repre- 
sentation and MF-DFA method argued as working better than other approaches [15]. Our main 
task is to investigate how multifractal properties of stocks change with variety of "surgeries" done 
on the initial real financial time series. We shall do this in following chapters. Our aim is to see 
what phenomena on the market influence its multifractal dynamics in the first place. In particular, 
we will focus on examining how multifractal picture of real time series is changed when one cuts off 
extreme events like crashes or rupture points, and how established trends or their absence influence 
the multifractal behavior of developed and developing markets. 

2 Multifractal noise in monofractal signal of finite length 

All statistical methods determining the fractal properties of time series give an exact result only for 

series of infinite length. Since this applies also to the MF-DFA scheme, one should know how the 
length of series is an important factor in discovering its multifractal structure. In fact, finite signals 
assumed by their construction to be monofractals, reveal some artificial multifractal structure. This 
is because big fiuctuations in finite time series, measured within MF-DFA, seem to be more rare 
than for longer or infinite time series. It leads to smaller h{q) value for g — > oo in finite series than 
for infinite ones, and finally, to some artificial multifractal structure of finite time series. This kind 
of multifractal noise should be subtracted first as a background infiuencing any real property of 
monofractal or multifractal time series of finite length. 

To explore this issue we must have a "test series" with known fractal characteristic. We will use 
monofractal series generated by the random midpoint displacement (RMD) algorithm [22]. The 
resulting h(q) functions calculated within MF-DFA for three series of length L = 2^^ generated with 
input Hurst exponent values H = 0.3, 0.5, 0.7 respectively are presented in Fig. 1. 
One may notice a small degree of multifractality (short span of h{q)) for all three series. The series 
with higher H show a bigger deviation from monofractality. The non-monotonic behaviour of h{q) 
is reflected by the "twist" on top of the singularity spectrum (see Fig. 2). This effect will be also 
observed in following chapters where real flnancial data are analyzed. 

One can also see that the maxima positions of the spectra reflect the assumed Hurst exponents 
values. 

It is essential to find the infiuence of series length on the width of multifractal spectrum. To 
observe this, one has to calculate the highest and lowest value of a, as well as the position of the 
f{a) peak. This was done on a sample of 20 series each for four lengths L — 2^, 2^°, 2-*^^, 2^^ and 
for different H values. The width of /(a) spectrum was rapidly changing with L and seemed to be 
weakly dependent on H. Table 1 presents results found for H = 0.5. 

They confirm the significant narrowing of Act spectrum width when the length L of time series is 
increasing. Since we will deal later on with series lengths around L ~ 2^^, we may assume that the 
finite size effects should be around SfinUei^c^) ~ 0.20. 
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Series 


Spectrum 


Peak 


length L 


width A« 


position 


2^ 


0.73 


0.56 


210 


0.28 


0.50 


212 


0.21 


0.50 


214 


0.22 


0.51 



Table 1: An example of dependence between the length L and the singularity spectrum width Act 
for artificially generated RMD series with assumed input value of Hurst exponent = 0.5. 

Before analyzing real data wc also need to observe the scaling range of the -F(s, q) ~ s^*-*^^ relation. 
For the artificial data, this range spans almost the entire available time window lengths (see Fig. 
3). However, for the entire set of real data (later sections will also discuss parts of these data), the 
scaling range becomes shorter and decreases with increasing |g| values. 

The corresponding effect for the S&P500 is visible in Fig. 4. It is even more significant for positive 
and larger q what is shown in Fig. 5 (compare with Fig. 4). To solve this problem a careful separate 
study of all F(s, q) ~ s^'^'^^ relations and manual selection of the scaling range is needed. It was done 
in analysis of all financial data presented below. 

3 Long-term data analysis 

The first and most obvious approach is to make the long-term data analysis, i.e., to calculate the 
multifractal structure of the stock market indices in the whole available time horizon. Our analysis 
focuses on two indices: 

• WIG (Warsaw Stock Exchange)) - from April 16*'* 1991 till October 10*'* 2008 

• S&P500 (NYSE & Nasdaq) - from December 30*'* 1927 till September 3^<^ 2008 

They are examples of developing (Poland) and developed (USA) markets respectively. 
All time series were created from daily closing values of each index. The S&P 500 data were down- 
loaded from the financial web-site of Yahoo.com [18] and WIG data from a stock market web-site of 
Wirtualna Polska [19]. 

The singularity spectra /(a) and h{q) plots of all discussed time series (as well as their shuffles) 
are presented in Figs 6, 7. Surprisingly, all markets show unexpected non-monotonic h{q) behavior. 
This phenomenon, in addition to the problems with short scaling ranges (mentioned in the previous 
section), strongly suggests a disturbance in the multifractal structure of the series. Spectra of all 
three shuffled data sets are positioned around the expected location a = 0.5. These sets of data also 
show a monotonic behavior of the h{q) function, contrary to originally ordered data. The spectra 
width for original as well for shuffled data arc much wider than those of the corresponding RMD 
spectra and hence indicate the multifractal content in dynamics of long-term market data. Simul- 
taneously, the non-monotonic h{q) behavior seems to indicate non-stationary effects in these data. 
The latter case needs more detailed study partly provided below. 



4 



4 Mult ifr act al analysis of financial data in established trends 



It has been shown by a number of authors (see e.g. [5,20,21]) that positive and negative fluctuations 
in time series have different fractal properties. Oswi§cimka et. al. has shown in Ref.[5] that DAX 
data in increasing and decreasing trends has different multifractal properties for particular choice of 
periods Dec. 1997-Dcc. 1999 and May 2004-May 2006. It is quite natural to investigate whether these 
properties are general and if so, in what extent they are general for other markets. In particular, 
it is interesting to know how multifractal properties of established trends depend on the internal 
structure of these trends. We were interested in the beginning in positive (negative) long-lasting 
trends with no sub-trends in opposite direction. Wc tried to exactly repeat this way the analysis 
made for DAX in Ref. [5]. Such trends turned out to be very rare for S&P 500 and WIG indices. 
The selected fragments of both indices with required properties are shown in Fig. 8. They have both 
long-lasting positive trend (bullish phase) followed by a long negative trend (bearish phase) and no 
astonishing internal events. 

The resulting singularity spectra (see Fig. 9) show the expected shift of the bearish phase spectra 
towards higher values of a. Since these results were obtained from data with a much wider scaling 
range, they should be considered as much more reliable then previous ones presented for entire in- 
dices. It is also worth noting that the spectra are now much wider, which suggest a more complex 
multifractal structure. In the case of WIG this effect might be amplified by the short data set, 
nevertheless the spectra width of a RMD series with similar length is still narrower. 
This improvement of singularity spectra for well determined trends in both indices (i.e. quasi 
quadratic dependence /(«) from a) with respect to spectra obtained for the entire indices might 
suggest that the multifractality is strongly affected by non-stationarity of financial time series. To 
verify such hypothesis we checked what multifractal properties are hidden in subparts of entire in- 
dex. The S&P500 index is particularly good for this search because it is long enough comparing with 
WIG, thus containing much more amount of data. If one divides the entire S&P500 history into two 
regimes - one with small fluctuations (small volatility) and the second one with high fluctuations 
(high volatility) as shown in Fig. 10, a very intriguing result for singularity spectra is found in both 
regimes. Results of this search are plotted in Figs. 11, 12, where a comparison of singularity spectra 
and ^(g) dependence for both regimes with /(a) and ^(g) behavior for the entire index is provided. It 
is clear that both regimes show separately a monotonic /i(g) dependence leading to reversed parabolic 
shape of /(a) . It is not the case of the entire index. 

However, both spectra for two regimes are slightly different - the width of multifractal spectrum 
(and the width of hi^q) dependence) is larger for the first regime characterized by lower volatility 
and no opposite sub-trends inside. We checked the same happened for other, even shorter trends, 
whenever they are affected by well formed sub-trends of opposite direction or by extreme events like 
crashes, etc. This phenomena is shown on example of S&P500 time series without and with internal 
structure (see Figs. 13, 14). Fig. 13 presents an example of increasing trend with no visible distortion 
inside. We called it "good" data. Its dynamics is associated with nice multifractal properties - well 
established quadratic /(a) dependence and the monotonic /i(g) behaviour. Simultaneously, time se- 
ries of the same length, but collecting data from the very next period, has some extreme events inside 
and it contains opposite sub-trends connected with these events (see Fig. 14). We call them "bad" 
data. The scaling Fg(s) ~ s'**^^^ is worse here and the presence of "twist" in singularity spectrum 
/(a) confirms the non-monotonic character of /i(g) function and leads also to non-quadratic /(a) 
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dependence. The multifractal character of time series is here very much affected by the presence of 
abrupt events, i.e. crashes or data with respectively large fluctuations. 

5 Influence of extreme events on multifractal dynamics of 
stocks 

The sub-trends containing abrupt events can be artificially removed to see if such removal would 
improve the multiscaling properties of the new trend constructed after such "surgery". The applied 
procedure is simple. First, the time series with abrupt events is differentiated. Then, at the level 
of returns, these events are removed and the time series of remaining returns is integrated to get 
the new artificial evolution of financial index. Examples of such modification are visible for variety 
of events for both S&P 500 and WIG indices in Figs. 15, 16 and 17 respectively. It is remarkable, 
how much such surgery significantly improves multi-scaling properties of financial time series and 
restores the monotonic character of decreasing h{q) function. The width of singularity spectrum has 
also noticeably increased when abrupt events had been removed. 

6 Conclusions 

We presented a comparative analysis of multifractal properties of financial time scries built from 
stock indices of developing (WIG) and emergent (S&P500) markets using MF-DFA technique. We 
found that financial time series in a very long time horizon show in both cases very short multiscaling 
ranges and in turn, their multifractal image is somehow obscured. Therefore we analyzed in shorter 
time horizon the selected parts of financial data chosen due to their specific properties. 

The most important result is that division of time series into regimes of distinctly different (large 
or small) fluctuations around the main trend, as well as removal of abrupt events like crashes or sig- 
nificant opposite sub-trends, radically improves the multiscaling ranges and restores the monotonic 
behavior of h{q) dependence as well as the reversed parabolic shape of multifractal singularity spec- 
trum f{a). The width of singularity spectrum also increases noticeably in these cases. This leads us 
to conclusion that non-stationarity of financial time series significantly infiuences their multifractal 
character leading to behavior beyond our expectation for multifractal spectrum shape known from 
studies with artificially created mono- and multifractal signals. 

It is clear that further research of length effects as well as the effects of non-stationarity in real 
time series influencing multiscaling and multifractality is needed with hope for further practical ap- 
plications. 
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Figure 1: The h{q) plot calculated for 3 artificial RMD series simmulated with assumed input values 
H = 0.3, 0.5 and 0.7 respectively. 
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Figure 2: Singularity spectra of 3 artificial RMD series with assumed H = 0.3, 0.5 and 0.7 respectively. 
The peaks of the spectra reflect the pre-set H exponent in terms of the Holder exponent a. 



8 




5 5.5 6 6.5 7 7.5 




log(s) 



Figure 3: An example of F{s,q) ~ s'^^'^^ plot for a series with assumed H = 0.7 for q = —5 (top), 
q = (middle) and q = 5 (bottom) 





1 














1 











S&P 500 q=-5 


5 


5.5 


6 


6.5 


7 


7.5 


8 


1 1 1 1 




1 










S&P 500 q=0 


5 


5.5 


6 


6.5 


7 


7.5 


8 




1 














1 










S&P 500 q =5 


5 


5.5 


6 


6.5 


7 


7.5 


8 



Figure 4: Exemplary plots of the F{s, q) ~ s^^'^^ relation for the entire S&P 500 series with q = —5 
(top), q = (middle) and q = 5 (bottom). 
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Figure 7: The singularity spectra of the entire S&P500 closing day data and its shuffle. The bottom 
figure shows the corresponding h{q) plots 
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Figure 8: Chosen parts of WIG and S&P500 indices containing the bearish and bulhsh phases. The 
moment of the trend change, sphtting the series into two distinct phases, is marked as a vertical hne. 
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Figure 9: Singularity spectra of the bullish and bearish phase for the S&P500. The bottom figure 
shows respective plots for WIG index. 
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Figure 10: The entire S&P500 time series. The vertical hne denotes the border between two regimes. 
One can see two phases corresponding to low and high oscillations around the trend. 
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Figure 11: Singularity spectra of the first part of the S&P500 and for the entire signal. 
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Figure 12: Singularity spectra of the second part of the S&P500 and for the entire signal. 




Figure 13: The top figure represents an example of "good" data without abrupt events. The bottom 
left plot shows an example of the F{s, q) function scaling for this data. The bottom right figure 
shows the corresponding singularity spectrum. 
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Figure 14: The top figure represents an example of "bad" data taken immediately after tlie period 
shown in Fig. 13. On the bottom left an example of F{s, q) ~ s'^*^'^-' scaling for this data is presented. 
It shows much poorer scaling properties than one in Fig. 13. The bottom right figure shows the 
"twisted" singularity spectrum, a result of a non-monotonic h{q) dependence. 
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Figure 15: The marked region with vertical hnes on the left shows the extreme event that was 
removed from the original S&P500 signal. The singularity spectra of the original part of the S&;P500 
and the modified one is shown on the right. 




Figure 16: The crossed region shows the abrupt regime. On the right, the singularity spectra of the 
original WIG fragment and the modified spectra after the abrupt regime has been removed. 
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Figure 17: Another example of abrupt regime removed from WIG time series (left) and the corre- 
sponding change in singularity spectra (right). 
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